Interaction-induced current-reversals in driven lattices 
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We demonstrate that long-range interactions can cause, as time evolves, consecutive reversals 
of directed currents for dilute ensembles of particles in driven lattices. These current-reversals are 
based on a general mechanism which leads to an interaction-induced accumulation of particles in the 
regular regions of the underlying single-particle phase space and to a synchronized single-particle 
motion as well as an enhanced efficiency of Hamiltonian ratchets. 
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Introduction Ratchets involve nonequilibrium dy- 
namical processes which make it possible to realize the 
old desire of converting thermal fluctuations into directed 
currents (DCs) and therefore usable work. According to 
their generality and their relevance for devices like Brow- 
nian [H-Q and molecular motors [|[ but also for biophys- 
ical problems like the migration of bacteria Q or the 
cell mobility in cancer metastasis [j| , the need for an un- 
derstanding of the detailed properties and prospects of 
the mechanisms of ratchets has given rise to the develop- 
ment of a highly active research field. This ranges from 
the theoretical analysis of the generating mechanism of 

DCs 0-G3 

to the remarkable realizations of ratchets in 
setups as different as semiconductor nanostructur es II ill . 
Josephson junction arrays (l2j and optical lattices |l3l - ll5| 
(see also the review [Ty]). Recently, it has been demon- 
strated that lattices with a spatially dependent driving 
imply a tunable phase space [171 ] and enrich the physics of 
DCs with mechanisms allowing for the creation of travel- 
ing density waves (l8j and designable patterned particle 
deposition [lj|- All of the above investigations do not 
focus on particle interactions and indeed most works on 
interacting ratchets concentrate either on the stochas- 
tic or the overdamped deterministic case poU^ leav- 
ing a gap in the literature concerning the microscopic 
analysis of interacting (deterministic) ratchets. A topic 
of p articular interest are the so-called current-reversals 
[n[ E3, [H-[ll, i-e. the tunability of the orientation of 
DCs via system parameters or particle density. Very 
recently even a current-reversal occurring in the time- 
evolution was achieved by a time-dependent modulation 
of the asymmetry of the ratchet potential [29| . 
In the present work we demonstrate that long-range in- 
teractions in a dilute particle ensemble cause self-driven 
current-reversals without requiring time-dependent mod- 
ulations of external parameters. Wc analyze this surpris- 
ing phenomenon to be the expression of a more general 
mechanism based on the interplay of two-body collisions 



and the underlying driven single-particle dynamics. As a 
consequence, we obtain synchronized single-particle mo- 
tion as well as an increase of DCs. 

Setup We consider a system of N equally charged par- 
ticles in a one-dimensional lattice of laterally oscillating 
Gaussian potential barriers with amplitude a, frequency 
uj, height V and equilibrium distance L described by the 
Hamiltonian 
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Therein Xi,pi,m denote the position, momentum and 
mass of the i-th particle, fj(t) is the driving law of the 
barrier with equilibrium position j ■ L and a := q 2 /4ire 
represents the interaction coefficient, where q is the 
charge of the particles and e is the permittivity. Af- 
ter performing the scaling transformations (x, t, m) <-$■ 
(x 1 := x/a,t' := tui,m' — muj 2 a 2 /V), we may set 
a = u) = V = 1 without loss of generality. For the 
remaining parameters we choose m = /3 = 1; L = 10 
and simulate 125 ensembles each consisting of N = 80 
interacting particles by numerical integration of the cor- 
responding Hamiltonian equations of motion. Large scale 
numerical computations have been performed for the in- 
tegration of the resulting 160 coupled, nonlinear ordi- 
nary differential equations for t/T = 2 • 10 5 periods of 
the driving law. Even though our results are essentially 
independent of the initial conditions, our focus is on di- 
lute ensembles whose initially stored interaction energy 
and kinetic energy is small in comparison to V (for sim- 
plicity all particles are placed equidistantly with a spac- 
ing Lq := Xi+\ — Xj). The total energy concerning the 
interactions of the i-th particle with the others reads 



Ma(x») = (a/L )[Hi. 
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i] where H k = £(l/j) 

i=l 

is the fc-th harmonic number. Choosing Lq = 15L 
and a < 8.0 we have V ia ,(xi) < 0.06a < V. The ini- 
tial particle velocities arc randomly chosen in the low 
velocity region of the chaotic sea of the corresponding 
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single-particle phase space (SPPS) resulting from Eq. Q] 
for N = 1. For fj(t), we choose a harmonic driv- 
ing law and in order to break the symmetries neces- 
sary to obtain a DC (Q) phase shifts of period three 
(^u4>2,4>3) = (0, 27r/3, 47r/3); i+3 = </>,; are applied 
[17l |. It is important to note that the observed effects 
and the underlying mechanism are generally valid and 
therefore not restricted to a specific driving as long as 
the relevant symmetries are broken. 

Current-reversals and fast ratchets We now explore 
the impact of long-range interactions on the directed 
transport. Fig. Q] (a) - (d) show the time evolution of 
the mean position and velocity of the ensemble of par- 
ticles for three different interaction strengths which can 
be adjusted according to the above-mentioned scaling be- 
haviour of the system. In Fig.[T](a) which shows the long- 
time behaviour for a = we observe a directed particle 
current with a negative transport velocity vt ~ —0.28. 
In the presence of interactions, we first note, that the 
directed transport persists. Surprisingly, for the inter- 
action strength a = 0.8 the direction of transport is in- 
verted. For a — 8.0, the DC points again in the same 
direction as for the noninteracting case. Fig. [1] (b) shows 
that in the noninteracting case a constant transport ve- 
locity is reached already after a short time t ~ 10 2 . In 
contrast to this, for a = 0.8 the DC initially points in the 
same direction as for a = and possesses a similar veloc- 
ity for several hundred of periods of the driving law but 
for longer times the particle current continuously slows 
down and reverts its direction. Even more, we observe 
a magnitude of the reverted transport velocity which is 
larger than for the noninteracting case (Fig.[T](c)). Thus, 
we encounter two phenomena: An interaction induced 
current-reversal and an increase of the ratchet-efficiency 
('fast ratchet'). For a — 8.0 even two current-reversals 
can be observed (FigQ] (a,b)). The DC deviates earlier 
from its noninteracting counterpart than for a = 0.8, but 
inverts its direction for a second time. Most clearly, the 
situation is reflected in Fig. [1] (d): Due to their small 
initial velocities the particles remain for the first few col- 
lisions on their lattice sites and are shaken according to 
the force provided by the driving law, so that their mean 
velocity simply follows the latter. Then, for a = 0.8 and 
a = 8.0 the mean velocity increasingly deviates from its 
noninteracting counterpart. We remark that in spite of 
the significant computational effort leading to the results 
shown in Fig[T]the asymptotic behaviour for long times 
is not yet reached. Indeed, below we will argue why the 
asymptotic states are reached only for such long times. 

Particle accumulation mechanism We derive now the 
mechanism which is responsible for the current-reversals. 
The backbone of our analysis is the single-particle dy- 
namics in the absence of all interactions. Since our many- 
particle system is dilute, the impact of the Coulomb in- 
teraction on the dynamics represents, most of the time, 
a small perturbation of the interaction-free dynamics. 
However, as we shall see there are events, namely two 
body collisions, for which the Coulomb interaction plays 
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FIG. 1. Time-evolution of the ensemble average of the posi- 
tion (x)(t) (a), its magnification for short and intermediate 
times (b), the average velocity (v)(t) (c) also with semiloga- 
rithmic axis (d). Parameters: a = u = V = /3 = 1; L = 10; 
a = (blue/medium saturation); a = 0.8 (green/low satura- 
tion) and a = 8.0 (red/full saturation). 



a significant role and which are responsible for the ob- 
served phenomena. Let us therefore first analyze the 
structure of the SPPS via the stroboscopic single-particle 
Poincare surface of section (PSOS) shown in Fig. [5J The 
phase space is mixed, i.e. we encounter a large chaotic sea 
at low velocities with embedded elliptic islands and in- 
variant spanning curves confining the chaotic sea. The el- 
liptic islands correspond to synchronized particle-barrier 
motion in configuration space. Islands crossing the v = 
line belong to particles which are trapped between two of 
the Gaussian potential barriers, whereas the others be- 
long to particles ballistically moving in a given spatial 
direction. The first invariant spanning curve (FISC) lim- 
its the energy which is maximally achievable for chaotic 
trajectories. Motion on these invariant spanning curves 
involves energies which allow the ballistically flying parti- 
cles to traverse the barriers for any phase of their motion. 
We observe an apparent asymmetry w.r.t. v = in the 




FIG. 2. Stroboscopic Poincare surface of section of the single- 
particle phase space (N — 1) at times t/T £ N. Parameters: 
a = oj = V = P = 1; L = 10. Blue dots: chaotic sea. Red 
lines: regular orbits. 



PSOS (Fig. [5]), which reflects the breaking of the rele- 
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vant symmetries [8J and is responsible for the occurrence 
of directed transport in the noninteracting system. 
For the interacting many-particle system it is not pos- 
sible to visualize the underlying high-dimensional phase 
space and we focus on specific observables to analyze and 
understand the time-evolution. Let us first inspect the 
time-evolution of the velocity distributions. For a = 
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FIG. 3. Time-evolution of the velocity distribution for 125 
ensembles of 80 particles. The number of particles (color) is 
shown as a function of time (logarithmically) and velocity, (a) 
a = 0; (b) a = 0.8; (c) a = 8.0. Ensemble and parameters as 
in Fig. [[] 

(Fig. El (a)) we encounter a uniform distribution in regions 
where we have no relevant elliptic islands in the phase 
space (Fig. [5]) and dips at velocities where large islands 
are located. The sharp decrease at v ~ 3.0 and v ~ —2.7, 
for times t > 100 when the single-particle phase space is 
occupied uniformly, reflect the velocity-boarders of the 
FISCs. The occurrence of a directed transport with neg- 
ative velocity is predominantly due to the strongly re- 
duced density for velocities v > 1.9 corresponding to the 
large elliptic islands centered at v ~ 2.5 (Fig. (2). 
For a = 0.8 (Fig. [3] (b)) and for short times, the ve- 



locity distribution is similar to the corresponding dis- 
tribution for a = 0, reflecting the observed equality of 
the transport velocity in both cases. Subsequently, the 
fine structure of the SPPS, i.e. the dips in the velocity- 
distribution, dies out continuously. Simultaneously the 
particles start to accumulate at velocities v ~ 2 — 3 
and around v <~ —2.5. Since the accumulation of par- 
ticles around v ~ 2.5 is dominant for large times, we 
obtain a DC pointing in the opposite direction as for 
a = 0, i.e. a current-reversal is encountered. For 
a = 8.0 (Fig. [3] (c)) a faster broadening of the velocity- 
distribution compared to a = 0.8 is observed. From 
t ~ 200 on a particle accumulation develops in the veloc- 
ity region v ~ 2 — 3 (first current-reversal) and spreads 
out up to u ~ 3.9. Subsequently a second particle ac- 
cumulation sets in for v ~ (—3.5,-2.8) which becomes 
dominant (second current- reversal). The key observa- 
tion to understand the above-discussed effects is that 
the velocity regime (v ~ 2 — 3), for which the particle 
accumulation occurs in case of the first current-reversal 
matches the velocity regime of the large elliptic islands of 
the SPPS (Fig. [2]). The additional above-mentioned par- 
ticle accumulations occur for velocities which are directly 
above/below the FISC-vclocities of the SPPS. These ob- 
servations indicate, that the structure of the phase-space 
of the corresponding single-particle system is still of cru- 
cial importance for the many-particle dynamics. To ver- 
ify the particle accumulation in certain (regular) parts of 
the underlying SPPS we lay the PSOS (Fig© and the 
position-velocity-distribution of the ensemble for large 
times (t/T ~ 2 • 10 5 ) on top of each other. For a = 
(Fig. [4] (a)) all particles are situated in the chaotic sea 
of the SPPS and the regular parts are unoccupied. The 
same holds also for the short time behaviour (t < 10 2 T) 
for a = 0.8 (not shown). Considering the case a = 0.8 at 
t = 2 • 10 5 T we observe a vital accumulation of particles 
in both large elliptic islands centered at v « 2.5. This 
accumulation spreads out into the region of the invariant 
spanning curves of the SPPS. A further accumulation oc- 
curs directly below the lower FISC. Since the one in the 
large islands for positive velocities is dominant, the aver- 
age velocity has become positive. We therefore conclude, 
that both the current-reversal and the high velocity of 
the DC for a = 0.8 are caused by a particle accumu- 
lation in the regular islands of the corresponding SPPS 
which are inaccessible for chaotic trajectories in the non- 
interacting case. For a = 8.0, the particle accumulation 
in the large regular islands takes place for earlier times 
and it is responsible for the first reversal of the DC. Fig.fJ] 
(c) shows the corresponding superposition for long times 
(t/T ~ 2 • 10 5 ) and exhibits a dominant accumulation 
below the lower FISC of the SPPS which illuminates the 
second current-reversal with the DC pointing finally in 
negative direction. 

Let us explore the doorway process into the regu- 
lar parts of the SPPS and the subsequent trapping for 
long times. A typical island-entering event is shown for 
N = 2, a = 0.8 in Fig. [5j Before the two-particles en- 
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FIG. 4. Superposition of the single-particle Poincare surface 
of section (grey/green) and the distribution of particles in 
the phase space at t = 2 ■ 10 5 T (red/yellow). Parameters: 
a = u = V = p = 1;L = 10; a = (a); a = 0.8 (b) and 
a = 8.0 (c). 



counter a close collision (i ~ 3.57 • 10 4 see Fig. [5] (b)) 
their dynamics is chaotic. Thereafter one particle is in- 
side the regular island and the other particle exhibits 
chaotic motion. The regular motion persists until a fur- 
ther collision occurs: then the particle inside the island 
(red) can either penetrate deeper into the island or it 
might exit the island and proceed with chaotic motion 
again. The interaction-induced doorway to the regular 
parts of the SPPS is therefore inherently a two-particle 
process. Note that contact interactions instead of long- 
range interactions would only yield the exchange of the 
velocities of both particles and no island-entering events 
would take place. 

Discussion It is important to note that the localiza- 
tion onto the elliptic islands of the SPPS induces a syn- 
chronized ballistic particle-barrier motion which is stable 
with respect to small perturbations. This mechanism 
combined with the tunability of the regular structures of 
the SPPS and their possible change in time via exter- 
nal parameters [l?], EH, H3] provides the perspective of 
controlling the time-evolution of the direction and mag- 
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FIG. 5. The doorway process into a regular island: Extract 
of the time-evolution of the velocities (a) and positions (b) 
for N = 2 and a = 0.8 in the vicinity of a two-body collision, 
(c) Total rate of particle collisions (per time period T) for 
a — 0.8 (black) and t _1 ' 2 -fit (red/dashed) with an inset using 
a logarithmic timescale. 



nitude of DCs. 

Let us address the question of the balance of entrance 
and exit processes concerning the regular parts of the 
SPPS in order to illuminate the particle accumulation 
process. For a dilute system the mean interaction energy 
is small in comparison with the mean kinetic energy of 
the particles. As long as the particles are largely sep- 
arated in configuration space their Coulomb interaction 
provides only a weak perturbation to the single-particle 
dynamics resulting in general only in a weak perturba- 
tion of the regular motion. For the case of two particle 
collisions however, their strong interaction destroys the 
tori of the SPPS for a short transient time. A random 
perturbation could in principle change the single-particle 
dynamics with equal probability for into and out-of is- 
land processes. The crucial difference is, that a chaotic 
trajectory can approach in line with the single-particle 
dynamics the border of any elliptic island embedded in 
the chaotic sea, but particles confined to orbits inside el- 
liptic islands possess a minimal nonzero distance to this 
border in phase space. Consequently, small perturbations 
can be sufficient to bring a particle across the border in- 
side an elliptic island, but in order for a particle which 
is deep inside an elliptic island of the SPPS to leave it, 
a whole sequence of adjusted perturbations is necessary. 
For this reason, in the average more particles enter the 
regular parts of the SPPS than are leaving them. Wc 
have confirmed this by additional simulations for weaker 
interactions a = 0.08 for which also the smaller regular is- 
lands of the SPPS possess an enhanced occupation. Since 
the average distance between the repulsively interacting 
particles increases with increasing time, the collisional 
interaction-energy decreases and the doorway process to 
the islands becomes less efficient which is in line with the 
decrease of the acceleration of the DC observed in Fig. [1] 
(d). Fig. [5] (c) shows that the collisional rate increases 
for short times but thereafter decreases rapidly, reaching 
approximately one collisional event per 10 5 particles and 
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period T at t = 2 ■ 10 5 . For long times the decrease is 
well described by a i _12 -fit, i.e. it is slightly faster than 
linear. The latter leads to a stabilization of the parti- 
cle accumulation for very long time scales, probably also 
in the asymptotic regime. Let us finally briefly discuss 
the timescales of the dynamics. The ergodic filling of the 
SPPS (t < 10 2 , Fig.©, current-reversals (t ~ 10 3 -2-10 4 , 
Fig. [TJ and the decay of the collisional rate (t ~ 3 • 10 4 
Fig. 13(c)) happen at separate timescales. As the incom- 
ing particle current into a regular region of the SPPS 
notably increases with the surface of the latter (i.e. the 
corresponding particle density increases with the frac- 
tion surface/volume) we obtain particle accumulations 
first for the small islands, followed by the large islands 
and then beyond the FISC (Fig. 1 (b),(c)). This cur- 
rent increases with a and the population of the small 
and the large islands (Fig. [3] (b),(c)) as well as the first 



current-reversal happen earlier for a — 8.0 compared to 
a = 0.8. We note that the described mechanism applies 
to other long-range interactions and in particular even to 
short-range interactions, but with a significantly lower 
efficiency. Clearly, it is possible to design the current- 
reversals via the structure of the underlying SPPS and 
the interaction strength. 

The predicted current-reversals, enhanced ratchet effi- 
ciency and the induced particle synchronization should 
be observable in experiments for electrons in AC-voltagc 
or laser driven arrays of semiconductor heterostructures 
11 1 or cold ions trapped in time-driven lattice potentials 
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